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Translation invariance, topology, and protection of criticality 
in chains of interacting anyons 
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Using finite size scaling arguments, the critical properties of a chain of interacting anyons can be 
extracted from the low energy spectrum of a finite system. In Phys. Rev. Lett. 98, 160409 (2007), 
Feiguin et al. showed that an antiferromagnetic (AFM) chain of Fibonacci anyons on a torus is 
in the same universality class as the tricritical Ising model, and that criticality is protected by a 
topological symmetry. We now study rings of interacting anyons on the disc as well as the torus, 
and show that analysis on the disc necessarily yields an energy spectrum which is a subset of that 
which is obtained on the torus. For a critical Hamiltonian, this subset always corresponds to the 
scaling dimensions of local scaling operators which respect the topological symmetry of the system. 
Related considerations are also shown to apply for open chains. 

PACS numbers: 05.30.Pr, 73.43.Lp, 03.65.Vf 



I. INTRODUCTION 

The study of collective states of anyonic excitations is 
an exciting and yet relatively unexplored area of con- 
densed matter physics. The nontrivial exchange be- 
haviour of nonabelian anyons may be exploited for uni- 
versal quantum computationfiiSi with the simplest suit- 
able model being that of Fibonacci anyons. It has been 
suggested that, as the non-Abelian component of the 
k = 3 Zfc-parafermion Read-Reyazi state,'^ they may ap- 
pear in the fractional quantum Hall state with filling frac- 
tion 1/ = 12/5^ These systems are therefore presently of 
in tense theoreti cal and experimental interest. 

iFeiguin et al\ recently initiated the study of interacting 
non-Abelian anyons with the analysis of nearest neigh- 
bour interactions in Ref. 5, in which it was shown that a 
chain of Fibonacci anyons on the torus could be analysed 
by mapping it to an equivalent spin chain. This work 
was later extended to next-to-nearest neighbour interac- 
tions and SU{2)k anyons in Refs. [y and[7|. These papers 
identified numerous critical phases, and scaling dimen- 
sions of the local scaling operators were extracted using 
exact diagonalisation by shifting and rescaling the result- 
ing energy spectra of the Hamiltonians.'^'^ Local scaling 
operators are of interest as they may appear as perturba- 
tions of the critical Hamiltonian, and may be classified by 
whether or not they respect the topological symmetry of 
the system. For Fibonacci anyons undergoing an antifer- 
romagnetic (AFM) interaction, the authors of Ref . [5| show 
that this topological symmetry protects the critical sys- 
tem against all conformally relevant translation-invariant 
perturbations. 

In this paper we investigate the relationship between 
periodic chains of anyons on the torus and on the disc, 
and their mappings to equivalent spin chains. We show 
that while the natural definition of translation for a ring 
of anyons on the torus is closely related to that on the 
spin chain, the same cannot be said for the natural defini- 
tion of translation on the disc, and a Hamiltonian which 
is translation invariant on the disc will only be transla- 



tion invariant up to a defect on the corresponding spin 
chain. The spectra of the same local Hamiltonian acting 
on two periodic chains of anyons, one on a torus and one 
on a disc, will therefore not in general coincide. We fur- 
ther show that the energy spectrum obtained on the disc 
always constitutes a subset of the spectrum obtained on 
the torus, and that for a critical theory, the local scal- 
ing operators associated with this subset are precisely 
those operators which respect the topological symmetry 
defined in Ref. ^ We also show that similar consider- 
ations apply to open chains, where the spectrum of the 
theory, and for critical theories the inferred local scaling 
operator content, is once again affected by the topology 
of the surface on which the anyons are found. 



II. ANYONIC STATES AND OPERATORS 



Although many papers have been published which 
study the behaviour of anyonic systems on surfaces of 
various topologies r^"—'^"— little attention has been paid 
to how the diagrammatic formalism may be used to ex- 
plicitly develop the relationship between states on sur- 
faces of various genus. In this section we address this 
deficit, beginning by reviewing the origin and formulation 
of the diagrammatic representation of states and opera- 
tors for systems of anyons on surfaces of genus (sphere, 
finite disc, and infinite disc) in Sec. Ill Al This material 
may be familiar to many readers. However, we present 
it here in a manner intended to emphasise the relation- 
ship between anyon models and topological quantum field 
theories (TQFTs)r^ii^i2i"— as we will exploit this rela- 
tionship to generalise the formalism to surfaces of higher 
genus in Sec. IIIBI Wc will also explicitly examine the 
construction of the translation operator on surfaces of 
genus and 1, as this will prove important to the study 
of translation invariant local Hamiltonians on the disc 
and the torus in Sees IIIIIIIVI 



A. Anyons on surfaces of genus 
(disc, sphere) 

1. Diagrammatic representation of states 

A system of anyons may be considered to consist of a 
collection of localised quasiparticle excitations in a two- 
dimensional medium, for example the topological liquid 
of a Fractional Quantum Hall (FQH) state^^ii^^"^ In 
general, a system may be considered anyonic if its par- 
ticles may be described in terms of a Unitary Braided 
Tensor Category (UBTC). In this paper we will concern 
ourselves only with anyon models which may be defined 
on the torus, known as modular anyon models)^^— for 
which the properties of the quasiparticles admit descrip- 
tion in terms of both a Unitary Braided Modular Tensor 
Category (UBMTC) and a 2-hlD Topological Quantum 
Field Theory (TQFT)i^ii^i2i-^ of the Schwarz type-^i^s 
Each of the quasiparticle excitations, or anyons, may 
then be characterised by a label, or charge, which cor- 
responds to a label of the UBMTC. 

However, providing a full description of such a system 
is in general more complicated than simply cataloguing 
the value and location of each non-trivial charge. This is 
because specifying the individual charges of two anyons, 
a and h, does not necessarily uniquely determine the total 
charge of the pair (a x &). These total charges are con- 
strained by the fusion rules of the UBMTC, which may 
be written in terms of the multiplicity tensor iV^;, as 



(i) 



a, a, a 
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but when there exist nonzero entries in N^^^ such that 
multiple terms appear on the right-hand side of this equa- 
tion, the total charge of a and b may correspond to any 
of these values c such that A^^^ 7^ 0. To specify these 
products, we represent the state of a system of anyons 
by means of a fusion tree (Fig. [1]). Labels on the inte- 
rior edges of the fusion tree graph correspond to the total 
charge of multiple anyons. For example in Fig. [TJi), xi 
is the total charge of anyons ai and 02 together, X2 is 
the total charge of anyons ai, 02, and 03 together, and 
so on. The set of valid labellings of a single fusion tree 
constitutes an orthogonal basis for the Hilbert space of 
a system of fixed anyons, and a labelling is deemed valid 
if all fusion vertices correspond to processes associated 
with non-zero entries in the multiplicity tensor A^j^. In 
this paper we will normalise all fusion tree bases using 
the diagrammatic isotopy convention given in Refs. l36| 
and [33. For free anyons, the co-ordinates of each anyon 
must be specified in addition to the fusion tree. 

Although a single fusion tree does not explicitly state 
the outcome of all possible measurements, it is possible 
to convert between different fusion trees using procedures 
known as F moves and braiding (Fig. [5]). In construct- 
ing a fusion tree, we have imposed a (possibly arbitrary) 
linear ordering on the anyons of the system. An F move 
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FIG. 1. Some possible fusion trees for a chain of six anyons 
with charges a\ to ar, on the disc or sphere. Labels Xi de- 
note intermediate fusion products which may not be uniquely 
determined by the fusion rules, and labels m axe associated 
with vertices and serve to enumerate multiple copies of a given 
charge for anyon models having some A^^ > 1- Note that 
no vertex index is required for fusion to the vacuum state. 
Tree (ii) is constructed from tree (i) by means of an F move 
[Fig. m^i)], and tree (iii) is constructed from tree (i) by recog- 
nising, from Eq. [l] that fusion with the vacuum state I is 
trivial. Diagram (iv) specifies a state Jt/j) of n anyons on the 
disc or sphere. 



[Fig-Uti)] alters the structure of the fusion tree while pre- 
serving that linear ordering, permitting the computation 
of additional fusion products [e.g. X2 in Fig. [Ijii) is the 
combined charge of 03 and 04], while braiding [Fig. [5l^ii)] 
permits conversion between different linear orderings.— 
Using these two operations it is possible to determine 
the probability amplitudes of different outcomes when 
measuring the total charge of any group of anyons re- 
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a, a, a. 







(ii) 




FIG. 2. Manipulations capable of performing a change of 
basis on a fusion tree: (i) F move, (ii) Braiding. 





UBMTC 



TQFT 



FIG. 3. Anyonic quasiparticles carrying labels from the 
UBMTC (denoted x ) map to punctures in the manifold when 
the system is represented by a TQFT. 



gardless of the fusion tree structure on which the state is 
initially described. The tensors ("fi^aaa") 

and Ral"^'^ are specified by the UBMTC to which the sys- 
tem of anyons corresponds. 

While the associated UBMTC describes a system of 
anyons in terms of individual quasiparticles, an equiv- 
alent description may also be made in terms of the 
diffeomorphism-invariant fields of a Schwarz-type 2+lD 
TQFT. Here, the 2D manifold on which the anyons exist 
becomes the spatial manifold of the TQFT, with individ- 
ual anyons corresponding to punctures in this manifold 
(Fig. [11). The state of the TQFT may be specified in 
terms of the outcome of a complete set of commuting 
Wilson loop operators, whose expectation values may be 
identified with the labels of the UBMTC. In a TQFT, 
a pair of Wilson loop operators which are topologically 
equivalent necessarily constitute a measurement of the 
same observable. Furthermore, the outcome of a Wilson 
loop measurement which may be contracted to a point 
is necessarily trivial. Consequently, we may identify the 
expectation value of an appropriate Wilson loop opera- 
tor (of specified orientation, to allow for charges which 
are not self-dual) with measurement of the total charge 
on the anyons, or punctures, which it encloses. Where 
degeneracies exist (i.e. N^^^ > 1 for some a,b,c), the dif- 
ferent copies of a particular charge label in the UBMTC 
can be associated with different expectation values of the 
Wilson loop operator. We therefore see that we may map 
between the TQFT and the UBMTC fusion tree as fol- 
lows: 

First, perform a pairs-of-pants decomposition of the 
punctured 2D spatial manifold of the TQFT. Then, take 




FIG. 4. Construction of a fusion tree graph from a pair of 
pants. 



one specific pair of pants (or 3-punctured 2-sphere) , de- 
clare that this 2-sphere has an inside and an outside, and 
specify which is which. Extend this definition of inside 
and outside consistently over all pairs of pants (this is 
always possible for a non-self-intersecting orientable 2- 
manifold embedded in M.^). Having done this, construct 
the fusion tree by drawing lines inside each pair of pants 
as shown in Fig. |4l Now associate a Wilson loop operator 
which measures charge with each opening of each pair of 
pants, up to topological equivalence. Specifically, where 
two pairs of pants connect together, we find two charge 
measurement operators which are topologically equiva- 
lent and so only one of these need be retained. Each line 
of the fusion tree graph now passes through exactly one 
Wilson loop, and we label the lines of the graph with 
the outcome of these charge measurements. Finally, if 
there exist entries in the multiplicity tensor N^^ which 
are greater than 1, then it is also necessary to associate 
a degeneracy index with the fusion vertices to enumerate 
these outcomes, which are also specified by the outcomes 
of the Wilson loop measurements. 

So far, this fusion tree has been constructed in the 
space M'^ in which the 2D spatial manifold is embedded. 
When representing this three-dimensional construction 
on paper, it is customary to employ a diagrammatic con- 
vention whereby the 2D manifold on which the punctures 
exist is mapped onto a plane perpendicular to the page, 
and whose projection onto that page forms a horizontal 
line at the top of the fusion tree diagram (for systems of 
anyons on the sphere, this is achieved by mapping that 
sphere first to the Riemann sphere, and then to the infi- 
nite plane) . The vertical axis of fusion trees drawn in this 
way (e.g. Fig. [ij may then be interpreted as a possible 
history whereby the present physical state may be ob- 
tained from the vacuum (i.e. a state with no punctures) 
in the 2-l-lD TQFT, with the lines of the fusion tree cor- 
responding to world lines of the quasiparticles presently 
observed on the manifold. For this reason, a charge la- 
bel I is typically placed at the bottom of the fusion tree 
diagram, representing the initial vacuum state4i 

By adopting different pairs-of-pants decompositions of 
the spatial manifold of the TQFT, it is possible to re- 
cover all different fusion tree bases of the UBMTC. It 
is also possible to interchange the definitions of "inside" 
and "outside" when constructing the fusion tree from the 
pairs-of-pants decomposition, but for surfaces of genus 
this has no effect on the basis obtained. In Fig. [5] we 
give a simple example of the pairs-of-pants construction, 
showing the decomposition of a 6-punctured finite disc 
with trivial charge on the boundary which corresponds 




FIG. 5. A sample pairs-of-pants decomposition for a 6- 
punctured finite disc: the manifold is decomposed into pairs 
of pants by cutting along the dotted lines. When the charge 
associated with the boundary is I, the construction described 
in Sec. Ill ATI yields the fusion tree of Fig. \V[i). 



to the fusion tree of Fig. [IJii). 

We conclude this section with a couple of remarks 
about specific systems of genus 0: First, to extend the 
pairs-of-pants construction to surfaces having less than 
three punctures, such as the 2-punctured 2-sphere, we 
note that fusion with the identity label I is trivial. For 
such a system we may therefore freely introduce addi- 
tional trivial punctures to obtain a single 3-punctured 
2-sphere from which we construct the fusion tree. Simi- 
larly, lines carrying trivial charge may be freely added to 
or removed from any fusion tree diagram [e.g. to obtain 
Fig. [Ijiii) from Fig. [Iji)]. Second, we note that there 
exists an important relationship between the sphere and 
the finite disc. While the infinite disc is topologically 
equivalent to the Riemann sphere, the finite disc may be 
treated as the Riemann sphere with a puncture at infin- 
ity. The edges of this puncture then constitute the edges 
of the disc. The charge associated with the edge of the 
disc is measured by a Wilson loop of the usual orien- 
tation enclosing this puncture on the Riemann sphere, 
or equivalently, one of reversed orientation enclosing all 
other punctures on the disc. When the charge associ- 
ated with this puncture on the Riemann sphere is (and 
remains) trivial, we may delete the associated line from 
the fusion tree diagram, and therefore ignore the exis- 
tence of the boundary when studying anyon behaviour 
on the finite disc. Third, we note that in the study of 
lattice models with n sites, we may treat the system as 
always containing n anyons at fixed locations, even if 
some of these anyons have trivial charge. The states of 
these systems can therefore always be represented by a 
fusion tree with n leaves. The enumeration of the leaves 
of the fusion tree then corresponds to an enumeration of 
the lattice sites, and consequently for such a system it is 
not necessary to separately state the co-ordinates of the 
individual anyons. 







FIG. 6. Elimination of loops during evaluation of the inner 
product. The numerical factor given is appropriate to the 
diagrammatic isotopy convention. 



a fusion tree to obtain a splitting tree, taking the com- 
plex conjugate of all fusion tree coefficients c''i-'^2n^ g^^^j 
reversing the direction of all arrows on the tree. In this 
paper, we will prefer lower indices for the coefficients of 
a splitting tree, e.g. c'^ ^ . The inner product of two 
diagrams is then performed by connecting the leaves of 
the fusion and the splitting tree, subject to the require- 
ment that leaves which are connected represent anyons 
(or punctures) at the same location on the manifold, and 
that the charges of the connected leaves coincide. Where 
these conditions do not hold, the inner product of two 
diagrams is zero. Recall that the fusion tree is a 2+1- 
dimensional structure projected onto a two-dimensional 
page, and thus when performing this connection, both 
trees must be represented in equivalent projections. Con- 
version between projections may be achieved by a se- 
quence of appropriately oriented braids. 

Assuming that the inner product has not yet been 
found to be zero, then once the trees have been con- 
nected, F moves are performed, loops are eliminated ac- 
cording to the rule given in Fig. [51 and trivial punctures 
are removed, until the resulting diagram has been re- 
duced to a number. This number is then the value of 
the inner product. Extension to states represented by 
a weighted sum over multiple labelled diagrams follows 
from bilinearity. 



3. Diagrammatic representation of operators 

Now that we have presented the diagrammatic formu- 
lation for anyonic states and for the inner product, we 
are in a position to construct anyonic operators. Where 
these operators act on the entire system, the construc- 
tion is trivial as an operator is constructed in the usual 
manner, as a sum over bras and kets: 



d = J2o^j\A){i^j\ 



(2) 



2. Inner product 

Next, we introduce the diagrammatic representation of 
the dual space and the inner product. In the diagram- 
matic representation of the space of states, the conjuga- 
tion operation ^ is implemented by vertically reflecting 



For anyons the bra is replaced by a splitting tree, the ket 
by a fusion tree, and the coefficient bears indices corre- 
sponding to all labels on the splitting and fusion trees 
[e.g. Fig. [TKi)]. However, we may also wish to define op- 
erators which act only on a finite subregion of the disc. 
In the same way that the fusion tree specifies how all the 
anyons in the system may be obtained starting from the 
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FIG. 8. Application of the translation operator T^ to a state 
of n anyons on the disc. The untranslated state \tp) is given 
in Fig. lUiv). 
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(iii) 



(iv) 



FIG. 7. Examples of anyonic operators on the disc with n 
punctures, (i) A global operator, acting on the total Hilbert 
space of the system, (ii) A local operator, acting on r adjacent 
anyons. (iii) The braid operator, (iv) The periodic translation 
operator T^ for a ring of anyons on fixed lattice sites on the 
disc, closing away from the observer. 



vacuum state, in an appropriate basis we may interpret a 
portion of the fusion tree as specifying how all the anyons 
within a physically localised subregion may be obtained 
from a single initial charge. For example, in Fig. [iji), we 
see that charges ai, 02 and 03 are obtained by splitting 
an initial charge of X2, and in diagram (ii), charges 03 
and 04 are obtained from X2- We require that our op- 



erators respect superselection rules associated with the 
charge labels of the UMBTC, and consequently they can- 
not change this total charge, but their action within this 
region is otherwise unconstrained. A completely general 
local operator acting on r adjacent sites on the disc may 
therefore be written in the form of Fig. Eljii). As with 
the state of a system, the choice of fusion and splitting 
trees employed in this figure merely represent a choice of 
basis in which to represent the operator, and any alterna- 
tive choice would have been equally valid. We also note 
that this construction includes the definition of a global 
operator on the disc, as the special case r = n. 

Finally, we note that while any operator on the disc 
may be represented in the form of Fig. Eljii), it may fre- 
quently be advantageous to represent certain special op- 
erators in other forms. Thus, for example, while the braid 
operator corresponding to the oriented exchange of a pair 
of anyons may be represented in the form of Fig. El^ii) for 
r = 2, it is usually more convenient to represent it in the 
form of diagram (iii), from which its unitarity is obvious 
by diagrammatic isotopy. Similarly, consider a ring of 
anyons occupying fixed lattice sites on the disc. Exploit- 
ing topological invariance, we may construct our fusion 
tree such that these lattice sites lie in a line at the top 
of the diagram, and the closure of the ring is implicit, 
being either towards or away from the observer. If, for 
definiteness, we assume that the ring closes away from 
the observer, then we may expediently represent the op- 
erator corresponding to periodic translation by one site 
using the diagram of Fig. [H^iv). Note that this opera- 
tor may be constructed by composing a series of braids 
[Fig.[71[iii)], and also that it respects the interpretation of 
the vertical axis as a fictional timeline for the creation of 
the state, as the motions of the anyons under the action 
of this operator are strictly monotonic in time. When 
this operator is applied to a state, the resulting diagram 
then describes a process whereby particles are created, 
migrate to their initial lattice sites, and then all move 
one site periodically around the lattice (Fig. [S]). 




FIG. 9. Construction of the torus from the sphere by in- 
troducing two punctures, deforming the resuhing punctured 
sphere by migrating the punctures towards the centre, and 
suturing. 
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FIG. 10. Operator Pt (in the diagrammatic isotopy conven- 
tion). 



B. Anyons on surfaces of higher genus (e.g. torus) 

Having reviewed the diagrammatic formulation for sys- 
tems of anyons on the disc or sphere, we will now extend 
this formulation to surfaces of higher genus, by exploit- 
ing the association between modular anyon models and 
2-l-lD TQFTs. 



1. Diagrammatic representation of states 




FIG. 11. Alternative procedure to construct a torus from the 
2-punctured sphere. The procedure presented here is compat- 
ible with a fusion tree constructed "inside" the torus, whereas 
the procedure presented in Fig. [9] is compatible with a fusion 
tree constructed "outside" the torus. 



Extension to surfaces of higher genus is achieved by 
means of manifold surgery, performed on the punctured 
manifold inhabited by the 2-f-lD TQFT. We will be par- 
ticularly interested in a specific example, the n-punctured 
torus, but the techniques which we will develop are en- 
tirely general and thus may be applied to construct dia- 
grammatic representations for states of anyonic systems 
on surfaces of arbitrary genus. 

We begin by noting that the torus may be constructed 
from the sphere by introducing punctures at the north 
and south poles, then distorting the sphere so that the 
puncture at the north pole descends vertically, and the 
puncture at the south pole rises vertically. When these 
punctures come into contact, they are sutured (Fig. IH]). 

Now, we wish to repeat this process for a manifold on 
which there exists a TQFT. We recognise that through 
the use of Wilson loop operators, charge labels may be as- 
sociated with the punctures cn and as at the north and 
south poles of the sphere respectively. On the sphere, 
prior to performing the suturing, these observables are 
independent. On the torus, after suturing, they are 
topologically equivalent up to a reversal in orientation. 
Importantly, these observables may be computed purely 
from the fields on the path of the loop itself, and thus 
their calculation proceeds identically whether or not the 
punctures are sutured. From this we infer two important 
results. First, suturing of these punctures only yields a 
consistent TQFT on the torus if the values of all Wilson 
loop observables on the north puncture are the duals of 
the same observables evaluated on the south puncture. 
Second, the space of states for the TQFT on the torus 
is isomorphic to the space of states on the 2-punctured 
sphere subject to this constraint. In Fig. [10] we see the 
operator Pt which projects from the Hilbert space of the 
n-|-2-punctured sphere to a reduced Hilbert space isomor- 
phic to the Hilbert space of the n-punctured torus. If we 



now describe the Hilbert space on the n+2-punctured 
sphere in terms of a fusion tree in the region of M3 collo- 
quially described as "outside" the sphere (i.e. extending 
from the surface of the unpunctured sphere to infinity), 
then we may use the surgical procedure described to con- 
struct a fusion tree for the n-punctured torus. Bringing 
together and suturing the punctures at the north and 
south poles corresponds to bringing together the equiva- 
lent branches of the fusion tree to form a loop. 

It is important to note that when constructing the 
torus from the sphere by means of the surgery procedure 
described, one necessarily obtains a fusion tree in the re- 
gion which is again "outside" the torus. This is because 
the branches of the fusion tree on the sphere which ter- 
minate in the north pole and south pole punctures must 
close to form a non-trivial cycle around the torus, and 
this can only occur if the fusion tree on the sphere inhab- 
its the "outside" space. A fusion tree "inside" the torus 
may be obtained by the alternative procedure of first con- 
structing a fusion tree "inside" the sphere, lengthening 
the sphere into an open cylinder with the polar punctures 
at its ends, and then bending this cylinder around into a 
loop and suturing (Fig. fTTj) . 

Given the existence of this relationship between anyon 
models on surfaces of higher genus and anyon models on 
the sphere, we see that on surfaces of higher genus we 
may employ the pairs-of-pants decomposition approach 
to construct a fusion tree basis in precisely the same way 
as we did for the sphere. Some example fusion trees 
for the n- punctured torus are shown in Fig. I12f i)-(ii). 
The fusion tree for the unpunctured torus is given in 
Fig. fT^ hi). and may be obtained using the pairs-of-pants 
approach by introducing a trivial puncture on the torus, 
constructing the fusion tree (where the puncture with the 
inward arrow in Fig. 2] is sutured to one of the punctures 
with an outward arrow), and then deleting the line car- 
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FIG. 12. COLOR ONLINE, (i)-(ii) Fusion trees for systems 
of n anyons on the the torus. The corresponding bases are 
related by means of a series of F moves, and will be denoted 
Bi and B2 respectively. In fusion tree diagrams on the torus, 
we will often place a * inside loops which correspond to a 
non-trivial cycle on the torus, to remind that these loops con- 
stitute an important part of the description of the state and 
cannot be eliminated using Fig.[Sl That is, the * signifies the 
presence of a topological obstruction. Note that for a given 
state, the value of x„ is unaffected by changing between bases 
El and B2 and this is reflected in the labelling of the diagram, 
(iii) Fusion tree for the unpunctured torus, (iv) Measurements 
Wa and Wb are associated with non-trivial cycles on the torus. 



rying charge I which is associated with the trivial punc- 
ture. Again in these examples, each line on the fusion 
tree diagram may be associated with a particular Wilson 
loop operator in the TQFT. However, on this occasion in 
addition to the measurements associated with punctures 
on the surface of the torus, there are also two measure- 
ments associated with the non-trivial cycles of the torus 
[Fig. HHiv)], and in a given basis, only one of these will 
encircle a line of the fusion tree. For example, consider 
the torus with no punctures. The fusion tree may be 



constructed either "outside" the torus (in the region of 
M'^ which extends to infinity), or "inside" the torus (in 
the region of K3 which does not extend to infinity) . The 
labellings of these two fusion trees both constitute a basis 
of states, and they are related by means of the topological 
S matrix. 
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(Note that an anyon model can consequently only be con- 
sistently defined on the torus iff the topological S matrix 
is unitary. This property is the defining characteristic of 
a modular anyon model.) 

In one of these bases, the fusion tree is encircled by 
Wilson loop operator Wa of Fig. [T^ iv). and in the other 
basis, by operator Wb- Thus by describing a state in one 
of these bases, we specify the probability amplitudes for 
the outcomes of measurements around both non-trivial 
cycles of the torus. 

We note that for the torus without punctures, the fu- 
sion tree admits as many different labellings as there are 
species of anyons in the model. When the fusion rules 
of the UBMTC correspond to a group Q, the number 
of labels corresponds to the number of elements in the 
group, \Q\, and we shall employ this notation even when 
the fusion rules do not form a group. The Hilbert space 
of the unpunctured torus is thus |5|-dimensional. As an 
example consider Kitaev's toric code, which is commonly 
understood as exhibiting independent electric and mag- 
netic charges, and has fusion rules corresponding to the 
group Z2 ®1i2- In the present terminology, we identify 
each element of Z2 Cg) Z2 as a separate charge, which we 
will denote I, e, m, and em. In the language of electric 
and magnetic charges, I is the uncharged vacuum state, 
e corresponds to the presence of an electric charge, m to 
a magnetic charge, and em corresponds to the presence 
of both. For the toric code, all states without punctures 
are ground states, and thus on the torus the ground state 
subspace has dimension 4, or equivalently we may say 
that the ground state on the torus is 4-fold degenerate. 
Similarly the dimension of the Hilbert space of states on 
an unpunctured manifold of genus g can easily be seen 
to be |5|^, reproducing the well-known ground state de- 
generacy of 4^ for the toric code on a surface of genus 

9- 

Finally, we draw attention to charge Xn-i in Fig. ll2f ii). 
Due to the presence of the topological obstruction de- 
noted by *, the loop in this fusion tree is not subject to 



the delta-function constraints of Fig. [13] which prohibit 
the existence of tadpole diagrams on the disc. On the 
torus, charge Xn-i is not constrained to be I. 



2. Inner product 

We now introduce a process for computing the inner 
product on the torus, which is derived from the inner 
product on the sphere by means of the process of man- 
ifold surgery described in Sec. IIIB II This construc- 
tion generalises immediately to all orientable non-self- 
intersecting surfaces of higher genus. 

Consider the inner product {4'''^\ip'^) between two 
states lip"^) and IV''"^) on the torus. For each state in 
turn we reverse the construction given in Sec. lIIBTl cut- 
ting the torus so that it is transformed into a surface 
isomorphic to the sphere with punctures at north and 
south poles, and then mapping each state |V''^), IV"'"^) on 
the torus to an equivalent state IV'^}, IV''^) lying within 
the support of Pj on the disc. As a notation convention, 
superscripts of T and D in this paper will be used to indi- 
cate that a particular state or operator lives on the torus 
or sphere/disc respectively. In contrast the t on Pr is 
written in subscript, and so is just part of the name we 
have chosen for this operator and does not denote the 
topology of the manifold on which the operator exists. 
The inner product between two states on the torus is 
now simply taken to be the inner product between the 
two equivalent states on the disc. 



(i) 



(vy^iV'T) = (v'°iv°). 



(6) 



We may summarise the computation of the inner prod- 
uct on the torus as follows: First, the fusion and splitting 
trees are connected at their leaves, as described for the 
sphere, and any mismatch between charges results in an 
inner product of zero. If the inner product has not yet 
been found to be zero, then F moves and Fig. [HI are ap- 
plied repeatedly until the diagram is reduced to a sum 
of terms having the form shown in Fig. 1131 These are 
then evaluated as shown, to obtain the value of the inner 
product. 

It is instructive to compare this formulation of the 
inner product with that presented in Appendix A of 
Ref. l2Ci The form ulation of the inner product intro- 
duced by iKonig an d Bilgin similarly guarantees that the 
physically permissible unique labellings of the fusion tree 
of the punctured torus yield an orthogonal basis for the 
Hilbert space, and differs only in the normalisation fac- 
tors which must be associated with some of the diagrams 
(see Table H. 




(ii) 




<a: Wn 




(iii) 



■ ^aia[^xix'^ V «i X 



FIG. 13. Evaluation of the inner product on the singly- 
punctured torus, in the diagrammatic isotopy convention. 
(i) Opposition of torus fusion and splitting trees, (ii) Equiva- 
lent diagram on the sphere, (iii) Numerical value. Note that 
in the diagrammatic isotopy convention, mapping states on 
the torus to states on the sphere effectively amounts to the 
removal of two vertices, along with their associated numerical 

factors. This introduces a factor of ^/di^dj.' in step (ii). 



genus. We will begin with a general discussion, and once 
again will examine explicit examples on the torus. 

On a surface of genus g, operators may correspond to 
physical processes acting either on the entire manifold, 
or on a finite subregion of the manifold. Where oper- 
ators act on the entire manifold, a completely general 
construction may once again be achieved by replacing 
the bras and kets of Eq. [5] with fusion and splitting tree 
diagrams for states of the appropriate genus. However, 
for operators acting on a finite subregion of the manifold 
the situation may be simplified somewhat. Momentar- 
ily neglecting the existence of punctures on the physical 
manifold, we examine the topology of the area of support 
of the operator. If this area of support now lies entirely 
within a submanifold of genus g' < g, then the opera- 
tor may be represented using fusion and splitting trees of 
genus g' . For example, if we consider an operator on the 
torus whose support lies within a region which is (again 
momentarily ignoring any anyons within it) topologically 
the unpunctured disc, then by locality we need only con- 
sider the portion of the fusion tree corresponding to any 
anyons which do lie within that disc. We then choose a 
basis where this portion of the tree connects to the rest 
of the fusion tree via only a single line, such that if we 
excised this disc from the manifold as a whole, that line 
would describe the charge on the boundary of the disc. 



TABLE I. Inner products of unnormalised diagrams on the 
1-punctured torus, for Fibonacci anyon statistics. Labels ai, 
xi, and x'l refer to diagram (i) of Fig. 1131 and a[ is set equal 
to ai. All inner products not listed below are zero in both 
conventions. 



3. Operators on surfaces of higher genus 

As with the sphere, we will now address the diagram- 
matic representation of operators on surfaces of higher 



xi, ai, x'l 


Convention of Sec. IIIB 2 


Convention of Ref. 20 


1,1,1 
r, I, T 

T, T, T 


1 


1 
1 



a, a, a, a. a, a. 



x,...x^ 

U,...LC 







\/ 



Xj...X^ 




FIG. 14. An operator O acts on a region of the 6-punctured 
torus which contains two anyons, and which, if these anyons 
were not present, would be topologically an unpunctured 
disc. The above diagrams represent an expression of the form 
lip'"^) = O^IV"^), and the basis on the torus has been chosen 
for convenience. 



We may now represent the operator in the form of an op- 
erator on the disc as described in Sec. IIIA3| and apply 
it by connecting it with the relevant portion of the fusion 
tree, as shown in the example of Fig. [141 We will, how- 
ever, make one warning: If the operator shown in Fig. 1141 
were to be applied to anyons ag and ai, then transfor- 
mation from the basis shown into one in which ag and oi 
were adjacent would require the application of the peri- 
odic translation operator, which is a non-local operator 
and will be discussed in Sec. Ill B 41 

Extension of this approach to surfaces and operators 
of higher genus is straightforward. 

There exists one further observation to be made with 
respect to surfaces of higher genus. Much as the torus 
admits operators of genus and genus 1, a surface of 
genus g will admit operators whose support is a region of 
genus g' , for any g' < g. We are not aware of any notation 
convention for the description of such operators, and on 
surfaces having genus higher than 1 , there is the potential 
for ambiguity as a given operator diagram of genus g' may 



refer to any of 




different physical processes. We 



suggest it may be appropriate to distinguish the different 
non-trivial cycles using different symbols in the manner 
of the * in Figs. [12] [HI but we do not develop such a 
system formally at this time, and note that ambiguity 



can always be avoided by writing such an operator in 
terms of full fusion and splitting trees on the surface of 
genus g. 

4- Periodic translation on the torus 

We now consider a specific example system which will 
be of interest in Sec. IIIII Suppose we have a system 
of anyons on a torus, arranged on a periodic lattice 
such that this lattice encircles either the large or the 
small non-trivial cycle of the torus 4^ As these two situ- 
ations are topologically equivalent, we choose it to encir- 
cle specifically the large non-trivial cycle with no loss of 
generality.*'^ How can we, in the diagrammatic notation, 
most efficiently represent the process of simultaneously 
translating each anyon around the torus by one site? 

If we construct our fusion tree "inside" the torus, ad- 
vance each anyon one site around the torus, and project 
onto the page, then the periodic translation operator is 
seen to act on the torus as shown in Fig. [TSj We may 
therefore write the torus translation operator simply as 



fT^ 



^ 



(7) 



Note the presence of the star in the diagrammatic rep- 
resentation of the translation operator, which indicates 
that the periodic translation passes around a non-trivial 
cycle on the torus, as opposed to merely cyclically per- 
muting the anyons locally on a disc-like region of the 
torus-shaped manifold. 

Implementation of cyclic permutation on the torus in 
the "inside" basis poses an interesting challenge. In con- 
trast with the cyclic translation operator on the disc [T^, 
Fig.[2Uiv)], the operator T'^ can not be constructed from 
local operations by composing a series of braids. Instead 
we must introduce another new operator, given in Fig. 1161 
which we will call the modified translation operator, T^. 
By diagrammatic isotopy, we see that the action of this 
operator T^ in basis Bi is to cyclically permute the de- 
grees of freedom oi, . . . , a„ and xi, . . . ,Xn- We may fur- 
ther use diagrammatic isotopy to redraw TjJ in the form 



fT 




(8) 



and this may be rewritten using Fig. [H^ii) as 



fT _ 

^M — 






^■ 



(9) 
Comparing with Eq. [71 we see that T^ and T^ differ 

only by the charge-dependent phase i?"""" . We therefore 
conclude that the application of the periodic translation 
operator on the torus to a labelled fusion tree in basis 
Bi constructed "inside" the torus is equivalent to multi- 
plication by R^"""" followed by cyclic permutation of all 
anyon indices ai, . . . , a„ and internal indices xi, . . . , a;„. 
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all indices , \-l/2 /, \-2 

{da,... da J (d,J 



(ii) 



Qj a.^ Oj 



ci„ > a„ , o„ 



a„ a, 0, a 





T 



T 



^ 



a„ a, 0, a 




FIG. 15. The translation operator T"""" on the torus, (i) represented as an operator in the bra-ket form of Eq. [2j and (ii) repre- 
sented as a mapping between a state \il)) and the translated state \ip') — i\ip)- The presence of the star in the braided portion 
of the diagram indicates that the periodic translation passes around the same non-trivial cycle on the torus as the loop labelled 



To construct a diagrammatic representation of the pe- 
riodic translation operator in an "outside" basis, we will 
proceed somewhat differently. This time, let us begin 
with a state written in basis B2 and constructed in the 
region "outside" the torus. First, we map this state to 
the n-|-2-punctured sphere. We then perform a further 
mapping of this sphere to the infinite plane, to obtain the 
situation depicted in Fig. WTi 'i) where the arrangement of 
punctures is shown on the plane of the page. Introducing 
a fusion tree for this arrangement of punctures, as shown 
in Fig. rTTT ii). it is easy to construct the appropriate trans- 
lation operator on the infinite plane [Fig. [TTT iii)]. This 
operator maps states in the basis of Fig. [ITfii) to states 
in the basis of Fig. fTTl iv). in which the anyon a„ is explic- 
itly braided around the south polar puncture. The equiv- 
alent operator on the torus in basis B2 is given in Fig. [THl 
where anyon a„ is seen to braid through the loop which 
carries the flux through the torus. This may also be in- 
tuitively understood by explicitly constructing the fusion 
tree in the "outside" space, as shown in Fig. [THl and ob- 
serving that during periodic translation of the punctures, 
one anyon is necessarily threaded through the loop of 
the fusion tree. Interestingly, and in contrast with bases 
constructed "inside" the torus, the translation operator 



for a basis "outside" the torus can be implemented en- 
tirely in terms of local operations once the state has been 
mapped to the equivalent sphere. This approach can not 
be applied to "inside" bases, as reversing the surgery pro- 
cess given in Fig. [Til involves cutting the ring of anyons, 
and this leaves the process of periodic translation on the 
equivalent sphere undefined. 

Once again, notice that in either basis, the translation 
operator on the torus respects the arrow of time of the 
associated 2-flD TQFT: In Figs. [H] and HH the trajec- 
tories of the punctures are all monotonic in the vertical 
direction. 



5. Topological symmetry operators 

Finally, we will find it useful to introduce one more 
class of operator on the torus which admits a special 
graphical representation. Consider now a torus with 
a ring of punctures around the large non-trivial cycle. 
These operators, which we will denote Y^ ^ describe a 
process whereby a pair of anyons carrying charges b and 
6 are created from the vacuum, travel around opposite 
sides of a non-trivial cycle on the torus coplanar with 
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(i) 



\/ 




FIG. 16. Definition of an operator 2^ on tfie torus, 
cyclically permutes the degrees of freedom ai,...,c 
xi, . . . ,x„ in basis Bi. 



which 
n and 



the ring of punctures and without braiding, and then an- 
nihilate back to the vacuum. Expressed as a map from 
a state |V') to a state \tp') where \tp) and \ip') are written 
in a fusion tree basis in the "inside" space, an operator 
Y^ may be written as shown in Fig. ISOT i). If we now 
re-express the state lip) in a basis constructed "outside" 
the torus using Eq. [5l operator Yj^ then takes the form 
shown in Fig. ISOT ii) . Using the identity 




(ii) 



''■^ 'i'. 1 o„ 




(iii) 



(iv) 



«„ a. a„ a, a, a 




db 



Sab 

Sib 



(10) 



we can see from Fig. ISOT ii) that Y^ will have eigenvalues 
Sba„/Sbi in the physical portion of the Hilbert space, and 
a state will be an eigenvector of ly iff it is not in a 
superposition over label Xn- 

In the present notation, the Y operator employed in 
Ref. [5| would be denoted Y!^ and is constructed in the 
"inside" basis, as per Fig. I^i). 



FIG. 17. (i) Aerial view of n -I- 2 punctures on the infinite 
plane equivalent to n punctures in a non-trivial ring on the 
torus. The fusion tree on the plane imposes a linearisation 
on these punctures, and we may choose this to be as given by 
the black line, (ii) The corresponding fusion tree. We may 
assume this fusion tree to inhabit the curved plane obtained 
by extending the black line of diagram (i) into the plane of 
the page. The grey arrows in (i) indicate the process of peri- 
odic translation of the anyons on the ring. Note that during 
the process of translation, one anyon crosses the plane of the 
fusion tree while passing between the punctures ajv and as- 
This is reflected in the periodic translation operator, labelled 
(iii). Application of the operator (iii) to states in the form of 
fusion tree (ii) yields states expressed in the fusion tree basis 
of diagram (iv) . 
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FIG. 18. Periodic translation operator around the larger non- 
trivial cycle of the torus, expressed in basis B2 constructed in 
the "outside" space, and represented as a mapping between a 
state \ijj) and the translated state l-^') = J |^). 



III. PERIODIC BOUNDARY CONDITIONS 

In this section we will consider periodic chains of 
anyons first on the torus (Sec. IIII C|) and then on the 
disc (Sec. IIII Dj) . We further specialise to models in which 
every site in the chain carries a fixed, identical charge, 
and consequently in this section, and also in the next, 
we set oi = a2 = ... = a„. For each topology (torus 
and disc) we will introduce a translation-invariant local 
Hamiltonian written as a sum of local terms, and take as 
a specific example the AFM nearest neighbour interac- 
tion for a chain of Fibonacci anyons. We assume that an 
operator which is local acts only on a disc-like subregion 
of the manifold (i.e. if it is an operator acting on the 
torus, it does not include any non-trivial cycles). Con- 
sequently such an operator may be written in terms of a 
fusion tree defined on the disc, as per Fig. El^ii). 

To express the states of our system, we shall use the 
basis given in Fig. [ijiv) on the disc, and in Fig. [T2lii) 
on the torus. The ring of punctures will be taken as 
encircling the large non-trivial cycle of the torus, and 
the fusion tree is constructed in the "inside" space. 

For the Fibonacci AFM interaction, which is a nearest 
neighbour interaction, all terms of the Hamiltonian take 




FIG. 19. In this diagram we see the process of periodic trans- 
lation represented schematically on the actual toroidal mani- 
fold, with the fusion tree visible in the "outside" region of R3. 
It is seen that, in bases constructed in the "outside" space, 
periodic translation on the torus threads an anyon through 
the non-trivial loop of the fusion tree. 



the form of Fig. [7l[ii) for r — 2. For clarity, from this 
point forwards we shall only provide explicit treatments 
for nearest neighbour Hamiltonians, though most of the 
arguments and techniques presented readily generalise to 
r > 2. 



A. Mapping to spins 

To study these systems, we will employ the technique 
described in Ref. [a, whereby the degrees of freedom for 
a one-dimensional system of anyons with fixed charges 
and nondegenerate fusion rules may be mapped to a spin 
chain. In the bases of Fig.[ljiv) and Fig. [121 1) (basis Bi), 
for a system of n anyons the Hilbert space of the system 
is spanned by the p free parameters of the fusion tree, 
xi . . . Xp, where p = n — 3 on the disc and p = n on the 
torus, and these p free parameters may be identified with 
a spin chain of local dimension \Q\. 

Because the Hilbert space of the spin chain is larger 
than the Hilbert space of the anyon chain, the Hilbert 
space of the spin chain is restricted to admit only those 
states which correspond to valid fusion trees under the 
anyonic fusion rules. 



B. Periodic translation on a chain of spins 

We note that there exists a special relationship be- 
tween the process of periodic translation on a chain of 
spins and periodic translation on a ring of anyons encir- 
cling a non-trivial cycle of the torus. Under the mapping 
of Sec. IIIIA| each degree of freedom xi, . . . , x„ on the 
torus is mapped to a site on the spin chain, and periodic 
translation on the system of spins, which we will denote 
T^, cyclically permutes these labels by one place. For a 
state satisfying ai = a2 ==... = a„, an equivalent effect 
may be obtained for the fusion diagram of the torus by 
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(i) 




(ii) 



Y,^ 



\l/ 





Y^- 



\l/ 



i*. 9 <3!„, a„ 




FIG. 20. Operator Y^ acts on a state on the torus. In (i), the anyon pair b and b travel around the large non-trivial cycle, and 
the fusion tree is constructed in the "inside" space. In (ii), the anyon pair still travel around the large non-trivial cycle, but 
the fusion tree is constructed in the "outside" space. 



applying the operator T^ discussed in Sec. Ill B 41 Fur- C. 
thermore, for fixed anyon charges oi, . . . , a„, the factor 
(i?°"°") in Eq.|9]is a constant. Consequently, for a ring 
of fixed, identical anyons on the torus, periodic transla- 
tion of the ring of anyons is equivalent up to a phase to 
periodic translation on the associated spin chain, and for 
Fibonacci anyons this phase is givenii by RC^ — e^'^'/^. 
For any operator O on such a chain, we therefore have 
the identity 



r'^Of^t ^fTQ^Tt ^fSQfS-\^ 



(11) 



Hamiltonian with periodic boundary conditions 
on the torus 



We now introduce a translation-invariant anyonic 
Hamiltonian, H^'^'^. The superscripts A, P, and T indi- 
cate that the Hamiltonian is anyonic, periodic, and con- 
structed on the torus respectively. For r = 2, we may 
write 



n-l 



n 



(12) 



Some care is still required in the computation of mo- 
menta of translation-covariant states, as when transla- 
tion by one site on the torus introduces a phase of e'^ , 
translation of the equivalent state by one site on the spin 
chain will introduce a phase of e''' =(7?""°") e'^ . In 
previous work^ this complication has been avoided by 
following a convention whereby momenta are computed 
relative to the ground state. In the present paper, how- 
ever, we will retain the full phase shift for didactic clarity. 



where local operator hf'i^i acts on lattice sites i and z-l-1, 
and takes the form of Fig. [^T] Unless otherwise stated, 
the evaluation of position indices such as i -I- 1 is assumed 
to be periodic in the range 1 . . .n, so (for example) site 
n -|- 1 is identified with site 1. 

As a specific example, we will consider the AFM in- 
teraction on the golden chain, for which all anyons a^ on 
the lattice are constrained to have charge t. Because the 
charges a^, a^+i, a,^, and a^_^^ are fixed and there are no 

degeneracy indices, we may denote the elements of hf-^j^^ 
by [hY'ij^i) where x corresponds to the fusion product of 
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, . , Til 

x.u.u' u 

A 

a. Oj'^j 

FIG. 21. Form of the two-site local operator used as a term in 
the local Hamiltonians H^''^''^ ([12} and _ff^'P'° 
ai, flj, di+i, and a^^i are assumed to be fixed. 



2T|) . Charges 



the Fibonacci anyons on sites i and z+l respectively. The 
AFM Hamiltonian favours the fusion path txt -^ 1 , and 
we therefore assign (hf-^ ^i)^ = — 1 and (hf^i+i) = 0. 

As demonstrated by iFeiguin et al\ ^ for the golden 
chain we may construct a three-body operator h on the 
spin chain whose action is locally equivalent to the two- 
body operator h^ on the system of anyons. To do so 
we introduce the spin chain equivalent of applying an F 
move at sites i and i -I- 1 of basis Bi, for i < n: 



F^-i.t.z+i\^^~i^^^^+i) 






(13) 



where indicates an operator acting on the spin chain. 
We then write 



?,S - { pS V f/S pS 



where 



/if|...,x.,...) = (/it,+i)..l 



..). 



(14) 



(15) 



For a chain of Fibonacci anyons, Oj = r for all values 
of i, and the above construction yields hf_i ^ ^_^_i corre- 
sponding to /ifj_|_j for all i < n. To construct the final 
term we will exercise some caution, as a„ and ai are 
presently located at opposite ends of the fusion diagram. 
We therefore begin with the expression 



/,A _ fTtA yiTt 

"■n,! — ^ ""n-l.n^ 



(16) 



Since we have fixed the charges of all punctures ai to be 
equivalently r, we may apply Eq. [TT]to obtain 

As T^ is equivalent to translation on the chain of spins, 
T^, we see that 



"n-l,n,l — ^ "'n-2,n-l,n^ 



(18) 



as might have been expected. The total spin chain Hamil- 
tonian may therefore be written 



H 



S,P,T _ V^ 

i=0 
n 



(fs)'(/.L,3UTSt 



(19) 






on a periodic spin chain of length n. 

For Fibonacci anyons with nearest neighbour inter- 
actions, Hamiltonians (fT2|) and (|19l) are both quantum 
critical Hamiltonians. The low energy properties of such 
systems are described by a conformal field theory (CFT), 
and the scaling dimensions of the local primary fields can 
be extracted from the low energy spectrum of the trans- 
lation invariant critical model on a finite lattice with pe- 
riodic boundary conditions 4i^ The results from exactly 
diagonalising ijS,p,T ^^^ AFM Fibonacci chains of lengths 
24 and 25 are presented in Table [Hi The energy eigen- 
values have been shifted and rescaled to give the scal- 
ing dimensions of the corresponding CFT, which for this 
Hamiltonian is the minimal model associated with tri- 
critical Ising model, M{4:, 3). 

For each scaling dimension. Table |ll] also gives a pa- 
rameter referred to as the "flux through the torus" . If we 
write our states on the torus in basis i?2 constructed in 
the "outside" region of M3, then this is simply the value of 
charge label a;„, which may be measured using the fam- 
ily of operators Y^^ . It is not difficult to see from Eq. [7] 
and Figs. [^ and BT] that any operator Y^ will commute 
with a Hamiltonian of local terms H^'^'^ on the torus, 
and thus for Fibonacci anyons we may simultaneously 
diagonalise ff^'P-T and Y^ . It therefore follows that we 
may associate every energy eigenstate with a correspond- 
ing eigenvalue of Y^ , which in turn corresponds to the 
measurement of a well-defined charge a;„ . 

For a translation-invariant Hamiltonian such as we 
have here, we may also assign a momentum to each state 
as shown in the dispersion diagram of Fig. [521 This di- 
agram clearly shows the distinction between (i) periodic 
translation on sites of the spin chain and (ii) periodic 
translation of anyons on the torus, with the difference 
between T^ and T^ resulting in a relative phase shift of 
i?j|"'^ = e^^^^^ as described in Sec. IIIIBI '^"^ Interestingly, 
the non-zero momentum of the ground state in Fig. 1221^ 11) 
implies that it is possible to distinguish between clock- 
wise and anti-clockwise translations of the ground state 
of the AFM golden chain. This observation is discussed 
in Appendix [XI 



D. Hamiltonian with periodic boundary conditions 
on the disc 



To construct a periodic translation operator on the 
disc, we recognise that the anyon sites in Fig. [ijiv) lie 
on a circle, which must be assumed to close either to- 
wards or away from the reader. Opting for the latter, 
we may define the periodic translation operator on the 
disc according to Fig. [71^ iv) . By inspection we see that 
translation may be implemented by means of repeated 
application of the braiding operator of Fig. [Tl^iii) , which 
we will denote Bf^^_^_^. Denoting the anyonic periodic 
translation operator on the disc of Fig. ^iv) by T^, we 
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have 



n-l 



j,D _ 



n^ 



A 



(20) 



TABLE II. Energy spectra for rings of (i) 24 and (ii) 25 Fi- 
bonacci anyons interacting via an AFM nearest neighbour 
interaction on the torus, shifted and rescaled to yield scaling 
dimensions for the associated conformal field theory. For even 
numbers of anyons, this is the minimal model M{4, 3), asso- 
ciated with the tricritical Ising model. For an odd number of 
anyons, we obtain the spectrum of M{4, 3) with a Z2 twist.— 
The scaling dimensions for odd numbers of anyons may be 
obtained from those for even numbers of anyons by fusing the 
corresponding scaling fields with the generator of this internal 
Z2 symmetry, the field e" , as noted in Ref . [3. 



(i) 24 anyons 


Numerics 


Prediction from GET 


Flux through torus 


0.0000 





I 


0.0750 


0.0750 (|j) 


r 


0.1989 


0.2000 (i) 


T 


0.8826 


0.8750 (|) 


I 


1.0622* 


1.0750 (|j + l) 


T 


1.1784* 


1.2000 (i + 1) 


T 


1.1841 


1.2000 (1) 


T 


1.8540* 


1.8750 (| + 1) 


I 


1.9469* 


2.0000 (0 + 2) 


I 


1.9843* 


2.0750 (|)+2) 


r 


2.0180* 
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FIG. 22. GOLOR ONLINE. Energy vs. phase diagram, where 
e' is the phase acquired by energy eigenstates on translation 
by one site, for the Hamiltonian of the golden chain with 
AFM interaction (i) made translation-invariant on the cor- 
responding periodic spin chain with p = 24, and (ii) natu- 
rally translation-invariant on a non-trivial ring of anyons on 
the torus with n = 24. Squares mark theoretical values for 
primary fields, and circles show selected descendants. Solid 
lines bound the energies from below. The system of spins 
yields a dispersion diagram identical to that of anyons on the 
torus, except that the momenta are shifted by an offset of 
-iln [{RVy^] = -¥ as described in Sec. HITBl 



We may now introduce a translation-invariant Hamilto- 
nian 



H 



A.P,D 



n — i 

n 



(21) 



* Eigenvalue is twofold degenerate 



where ^ denotes that this Hamiltonian acts on the disc. 
As before, we also introduce a spin chain whose sites 
correspond to the degrees of freedom of the anyonic fusion 
tree. However, on this occasion we use the fusion tree of 
Fig.[TJiv) as a basis of states, and the corresponding spin 
chain is of length p = n — 3. Near the centre of the chain 
we may construct the spin equivalent of an F move in the 
same way as for the torus P^ . but we see that operators 
near the end of the chain will act on a reduced number of 
sites. For example, an F move acting on sites a2 and 03 
acts only on spin variables xi and 2:2, mapping |a;ia;2) into 
Y.i^{FS^'' )xiS:i\S:iX2) ■ We will continue to denote the 
spin chain counterparts of these operators by Ff^^^ j_|_2, 
with the understanding that when evaluating Eq. [TJ] on 
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the disc, any indices Xq or Xp+i are to be replaced by 
charges ai and a^ respectively, and any indices X-i or 
Xp+2 are to be replaced by the vacuum charge I. We do 
not modify the spin chain, which continues to run from 
xi to Xp. This behaviour manifestly breaks translation 
invariance on the spin chain. 

For values of i sufficiently distant from 1 or n we may 
also map hf^ij^i onto a three-site spin operator as before, 

although this is now denoted ft.f_2 i-i i as it acts on spin 
sites Xi-2, Xi-i, and Xi. By using the extended defini- 
tion of F^i^i i_|_2 we may even write down spin operators 
equivalent to h^2^ ^2 3' ^n-2n-n ^^^ ^n-i.n- However, 
for h^ ^ we must introduce the spin chain equivalent of 

the anyonic periodic translation operator on the disc, T^ . 
To do this, we first construct the spin chain counter- 
part to the anyonic braiding operator given in Fig.[7l^iii). 
This is achieved by introducing a unitary operator Rf 
derived from the tensor R in Fig. EJii), which operator 
multiplies a state \xi) by a phase i?^^"''^''*^^. Using this 
we may write the spin chain equivalent of Fig. [Tjiii) as 






(^- -'' 



:,i+l,i+2J 



^i+l-'^^i,i+l,i+2- 



(22) 



As with F^ , the same special identifications for x_i, xq, 
Xp+i, and Xp+2 must be made when applying either R^ or 
B^ to a state. Using B^ we can define an operator T'^ on 
the spin chain which is equivalent to periodic translation 
on the lattice of anyons. 



P-t-i 

rp/S _ TT rS 
i=0 



(23) 



and thus compute the spin chain Hamiltonian which is 
equivalent to H^'^'^: 



p+2 



H^'"'°-E(^")XM(^"')' 



(24) 



The Hamiltonian ijS,p,D jg (,\(.Q^j-iy j^q^ translation invari- 
ant under the natural definition of translation on a pe- 
riodic spin chain. However, it does exhibit translation 
invariance under the action of the anyon-derived transla- 
tion superoperator T'S(.)j"St 

Away from the edges of the fusion tree, the action 
of T"S(.)j"St is equivalent to translation on the sys- 
tem of spins, such that (for example) r'^(/if 2 3)^'^^ = 
/i2 3 4- However, this does not hold where the transla- 
tion would yield an operator crossing between sites 1 and 
p. Instead, T'^ {hp_2.p~i,p)'^'^^ yields a two-site opera- 
tor acting on spin sites p — 1 and p, and it is neces- 
sary to apply T'S(-)T'St six times to map /ip_2,p_i,p into 

^i2 3i with none of the intermediate terms resembling 
a translation on the spin system of the original opera- 
tor hp_2p_ip- Nevertheless, the complete Hamiltonian 
satisfies T'S(ffS.P.D)rst = i^s,p,D^ 
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FIG. 23. COLOR ONLINE. Energy vs. phase diagram, where 
e'" is the phase acquired by energy eigenstates on translation 
by one site, for the golden chain with AFM interaction on the 
disc with n = 24. Squares mark theoretical values for pri- 
mary fields, and circles show selected descendants. Solid lines 
bound the energies from below. Note that phase shifts are 
in agreement with those computed on the torus [Fig. I22f ii')]. 
and are offset by ^ relative to the corresponding values on 
the spin chain [Fig. EW i)]. 



The results of exactly diagonalising JJS,p,d ^^^ given in 
Table IIII[ and a dispersion diagram is plotted in Fig. [23] 
for comparison with the torus [Fig. [2^ ii)]. As with the 
torus the energy spectrum has been shifted and rescaled 
to give the scaling dimensions of local scaling operators 
inX(4,3). 

This behaviour, where translation invariance exists rel- 
ative to an operator which is not the natural translation 
operator on the spin chain, has previously been observed 



TABLE III. Energy spectra for rings of (i) 24 and (ii) 25 
Fibonacci anyons interacting via an AFM nearest neighbour 
interaction on the disc, shifted and rescaled to yield scaling 
dimensions for the associated conformal field theory. For even 
numbers of anyons, this is the minimal model A1(4, 3), asso- 
ciated with the tricritical Ising model. For an odd number 
of anyons, we obtain operators from the spectrum of A1(4, 3) 
with a Z-i twist.— The scaling dimensions for odd numbers 
of anyons may be obtained from those for even numbers of 
anyons by fusing the corresponding scaling fields with the 
generator of this internal Z-i symmetry, the field e" , as noted 
in Ref . i5|. 
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* Eigenvalue is twofold degenerate 
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for certain 5't/(2)fc-invariant spin chain Hamiltonians by 
Grosse et al.^ It is now known that a relationship ex- 
ists between 5'C/(2)fc-invariant spin chains and chains of 
SU{2)k anyons, and although present research has con- 
centrated on anyo ns on the toru s, ^ it is nevertheless likely 
that the models of lGrosse et ali may similarly be mapped 
into interactions of SU{2)k anyons on the disc. The form 
of the anyonic translation operator also has practical im- 
plications for the restriction of the Hilbert space of the 
spin chain mentioned in Sec. IIII Al and these technical 
details are discussed in Appendix |B] 



E. Relationship between the torus and the disc 



Mapping between disc and torus states 



2. Conformal field theory with a defect 



By comparing Tables HI] and IIIII we see that on the 
disc, we compute scaling dimensions which correspond 
to those obtained for trivial flux through the torus. The 
reason for this may be seen by comparing Figs. [IJiv) and 
Fig. IT^ ii) . By restricting the flux through the torus a:„ 
in basis B2 to be I, we obtain a fusion tree identical to 
that of Fig. [TJiv) . The action of the translation-invariant 
Hamiltonian 



n — 1 



i=0 



(25) 



where X stands for T on the torus and D on the disc 
is therefore equivalent in both cases, and we obtain the 
observed correspondences between the energy spectrum 
of the torus and the disc. 

Away from criticality, the Hamiltonian is insensitive 
to non-local properties of the system. In the thermody- 
namic limit each energy level is therefore |t^|-fold degener- 
ate on the torus, with the degeneracy enumerated by the 
different values which may be taken by the flux through 
the torus. As the disc corresponds to the torus with the 
flux constrained to be I, the eigenvalues on the disc in the 
thermodynamic limit are the same, but nondegenerate. 

At criticality, the energy spectrum is no longer nec- 
essarily independent of the flux through the torus, and 
the degeneracy of energy levels on the torus may be bro- 
ken, with different flux sectors exhibiting different energy 
spectra. Nevertheless, the identification between the disc 
and the torus with trivial flux persists, and for critical lo- 
cal Hamiltonians applied both on the torus and on the 
disc, the disc necessarily exhibits the same energy spec- 
trum as obtained on the torus when the flux through the 
torus is constrained to be I. 

Recall now that each scaling dimension obtained from 
the energy spectrum is associated with a local scaling 
operator.Sj^ In Ref. [j, it is argued that the classification 
of scaling dimensions on the torus by the value of the 
flux through the torus translates into an equivalent clas- 
sification of local scaling operators, and that when in the 
ground state, only local scaling operators associated with 
a flux of I may cause local perturbations. 



In Sees . IIII Cl and HIl Dl we have seen that a spin chain 
of length p may be used, via appropriate mappings, to 
represent states of a system of either p anyons on the 
torus or p -I- 3 anyons on the disc. We also observed 
that each of these systems comes with its own definition 
of translation invariance. For the torus this corresponds 
(up to a state-dependent phase) to the natural definition 
of translation on a periodic spin chain, whereas for the 
disc this is given by operator T^ of Eq. (P5)l but never- 
theless corresponds (when applied to an operator using 
the adjoint action) to the natural definition of transla- 
tion invariance on the spin chain on sites sufficiently far 
from xi and Xp. It is therefore natural to interpret the 
difference between these two models as being equivalent 
to the introduction of a defect in translation invariance. 
We will show by construction that this defect has the un- 
usual property of being invertible. That is, there exists 
a second defect which, when introduced manually, will 
annihilate the original defect and restore the full spec- 
trum for a system of anyons on the torus, albeit a torus 
of length p — 3. 

We now construct a Hamiltonian on a disc of Fi- 
bonacci anyons which reproduces the spectrum of the 
AFM Hamiltonian on the torus. This Hamiltonian satis- 
fies translation invariance on the disc — i.e. it is invariant 
under the adjoint action of T° (PO)) — except for two local 
terms. These terms define a defect T). The spectrum of 
the resulting Hamiltonian for a ring of n anyons on the 
disc is equivalent to that of n — 3 anyons on the torus. 
If n is even, then 71 — 3 is odd, and as noted below Ta- 
ble UTJ the spectrum of a ring of an odd number of anyons 
on the torus is equivalent to that of an even number of 
anyons with a Z2 twist 4^ Thus the fusion of defect V 
with the Z2 twist constitutes the inverse of the defect 
which arises from application of the translation superop- 
erator f°(-)f'°'f. The Hamiltonian exhibiting defect V 
takes the form 



H 



A,P,D^T 



n — 3 



i+l 



n — 2,n— l,n.l,2 



''n-l,n,l,2,3 



(26) 
for a ring of n anyons on the disc, and details of its con- 
struction for the AFM or FM Fibonacci chain are given 
in Appendix [Cl 

As a consequence of the invertibility of the defect in 
translation, it follows that we may compute the spectrum 
of a system of anyons on the torus or the disc using a 
system of anyons also, in each case, either on the torus 
or on the disc. In this section we have shown how to 
extract the spectrum of a system of anyons on the torus 
from a system of anyons on the disc, by means of a local 
modification to the Hamiltonian. In Sec. IIII ET] we noted 
that one can extract the spectrum of a system of anyons 
on the disc from a system of anyons on the torus by 
means of a global operator restricting the flux through 
the torus to be I, and for completeness we now note that 
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FIG. 24. (i) Fusion tree used as a basis of states for the n-site 
Heisenberg spin chain. Note that in contrast to Fig. [Hiv), 
the total charge (or spin) is not constrained to be zero, and 
a tree with total charge Xn-i therefore represents a vector 
space of dimension dx^_j^- (ii) Diagrammatic representation 
of the Hamiltonian for the AFM spin-i Heisenberg model. 



this may be achieved by adding to the Hamiltonian an 
appropriate 1-site local operator acting on site x„ the 
spin chain. This 1-site operator applies an arbitrarily 
large energy penalty to states for which a;„ 7^ I. When 
this term is introduced, the resulting Hamiltonian has a 
spectrum equivalent to a system of n anyons on the disc. 
Understanding that systems of anyons on the disc and 
on the torus each admit their own natural notion of trans- 
lation invariance, and that these are equivalent up to a 
specific defect, provides another interpretation of the dif- 
ferent operator spectra of Tables [III and I II 1 1 The study of 
defects and boundaries in conformal field theory, both in 
the continuumS.!^ and on the lattice,— has shown that 
the presence of a local defect in a critical system may 
globally affect the scaling operator spectrum. In the 
present example, we have shown that the AFM golden 
chain on the disc is equivalent to the AFM golden chain 
on the torus with a defect, and that this defect eliminates 
the scaling operators whose scaling dimensions exhibit a 
flux of T in Table IHl and likewise eliminates their corre- 
sponding eigenvalues from the spectrum of the Hamilto- 
nian. 



F. Comparison with models on a spin chain 

In this section we have dealt with models such as the 
AFM golden chain, which are primarily defined on a ring 
of anyons. It is instructive to compare these with models 
such as the Heisenberg model, which is defined on a spin 
chain but which can also be expressed in the diagram- 
matic notation used in this paper. As an example, we 
will now consider the spin-^ AFM Heisenberg model with 
periodic boundary conditions, assumed to be constructed 
on a manifold which is topologically the disc. Because 
this model possesses SU(2) symmetry, we may represent 
it in the diagrammatic notation using a UBMTC based 
on the fusion rules of SU(2). States can be represented 
in the form of Fig. [MT i). and the Hamiltonian takes the 
form of Fig. [Mlii) ■ 



We may now analyse this system in two different ways. 
Either we may obtain the energy spectrum by exactly 
diagonalising the original spin chain, and compute mo- 
menta using the natural definition of translation on that 
spin chain, or we may write it in the diagrammatic nota- 
tion, and map this to a different spin chain as described 
in Sec. IIIIAI We would then compute the energy on the 
chain of fusion tree variables xi, . . . , Xn-i, and the mo- 
menta using translation operator T^ from Fig.[IJi). As is 
to be expected, the results obtained using these different 
methods agree to the limits of numerical precision. 

This observation has bearing upon the definition of the 
translation operator. If we use the definitions for T and 
r"^ given in Fig. [TJi), Eq. [71 and Fig. [THl then we find 
that the ground state of the AFM golden chain has non- 
trivial momentum. Suppose that, instead, we assume a 
momentum of zero for the ground state of the golden 
chain, and adopt 



T 



M — \-^l ) 



as the translation operator on the disc and 



(27) 






{Rr-) ^~y-y- -7 



^ 



® 

on the torus (noting that 3^ is the modified translation 
operator originally introduced in Sec. IIIB 41 correspond- 
ing to cycling of the fusion tree variables xi,...,Xn_i). 
For consistency we would then also have to use Eq. [271 
when working with the diagrammatic representation of 

the AFM Heisenberg chain, with i?j""^ = i?|^ = -1. 
However, the momenta obtained using this operator are 
inconsistent with results obtained by exactly diagonalis- 
ing the original spin chain, indicating that T^ ^ and not 
T]^, is the correct definition for the periodic translation 
operator on the disc. As we require that the torus with 
trivial flux be consistent with the disc, we also obtain 
that T"^, and not T'^, is the correct periodic translation 
operator on the torus. Thus study of the AFM Heisen- 
berg spin chain supports our claim that the ground state 
of the AFM golden chain has non-zero momentum, as 
observed in Fig. EWii) for the torus and Fig. [521 for the 
disc. 



IV. OPEN BOUNDARY CONDITIONS 

On 1-D systems with open boundary conditions the 
situation is somewhat simpler, but again some care must 
be taken as the fusion tree basis will again depend upon 
the topology of the quantum liquid. For example, as 
in Ref. y, one might choose to study the Hamiltonian 
corresponding to free boundary conditions on the torus. 



H 



A.F.T 



E^: 



A 



(28) 
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where F denotes free boundary conditions, which maps 
to the spin chain as 



zjS,F,T „ V^ ?s 



(29) 



(This may be contrasted with Eq. 1190 Similarly, one 
could place the same Hamiltonian on the disc: 



^ 'F,D _ y^/t^ .^^^ 



H 



n-1 
S,F,D _ Y^ f,S 



E^^ 



i — 2.i—l.i' 



(30) 

(31) 



i=l 



Once again the spectrum for the Hamiltonian on the 
disc is seen to be a subset of that on the torus ( Table ITV]) . 
and once again by means of appropriate modifications of 
the Hamiltonians, corresponding to alternative choices of 
boundary conditions, we may obtain either set of scaling 
dimensions on either topology. 



V. CONCLUSION 

In summary, this paper may be divided into two parts. 
In the first (Sec. [IT| we have drawn attention to the im- 



TABLE IV. Numerical results and CFT assignments for 
the smallest scaling dimensions on open chains of Fibonacci 
anyons of length n with AFM coupling and free boundary 
conditions: (i)-(ii) On the torus, Hamiltonian H°'^''^ 
(iii)-(iv) On the disc, Hamiltonian jf{°'^'^ 



portance of manifold topology in the study of anyonic 
systems. By appeal to the underlying TQFT, we have 
explained how to construct diagrammatic representations 
of anyonic states, operators, and the inner product, for 
surfaces of arbitrary genus, and we have done so explicitly 
for the torus, sphere, and disc. Applying these construc- 
tions to systems of fixed Fibonacci anyons (the golden 
chain), we have been able to construct explicit relation- 
ships between the energy spectra of equivalent models of 
interacting anyons on the torus, sphere, and disc. For 
critical systems, this implies an equivalent relationship 
between the scaling operators of the system. 

In the second part of the paper (Sees. IIIIHIVP we used 
these results to study the behaviour of an example sys- 
tem, consisting of a ring or chain of interacting Fibonacci 
anyons on either the torus or the disc. It has been shown 
that this chain is described by the same CFT as the tri- 
critical Ising model, and that on the torus its critical- 
ity is topologically protected.^ We have shown that that 
criticality is similarly protected on the disc. We further 
demonstrated that protection of criticality on the disc 
may be understood in terms of conformal field theory, 
where the system on the disc maps into a system on the 
torus with a defect in translation, and presence of that 
defect modifies the scaling operator spectrum which is 
exhibited (Sec.|IIIE2|. 

As a whole, this paper therefore presents the means 
to relate systems of interacting anyons on manifolds of 
differing topology, and applies this to examples using 
Fibonacci anyons. Insight is gained into the topologi- 
cal protection of criticality of these systems and into the 
robustness of this protection across surfaces of different 
genus, with equivalent protection shown to be exhibited 
on both the torus and the disc. 
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* Eigenvalue is twofold degenerate 
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Appendix A: Chiral symmetry and Fibonacci anyons 



In Sec. nil CI we noted that it was possible to distin- 
guish between clockwise and anti-clockwise translations 
of the ground state of the AFM golden chain, implying 
that the ground state is not reflection-invariant, i.e. that 
it is chiral. We will now explain this observation in more 
detail. 

There are four different Unitary Braided Modular Ten- 
sor Categories (UBMTCs) having Fibonacci anyon fusion 
rules, and they may be grouped into two pairs, the first 
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having i?^^ — e ^^^1'^ and R\'^ = e^'^'/^, and the second 
having i?[^ = e^""'/^ and i?^ = e-^'^'/s, (The members 
of each pair are distinguished by differences in signs in 
their F moves.) We will denote these UBMTC by Fib!^, 
where i £ {1,2} indicates the set of F moves being em- 
ployed, and ± refers to the sign on the exponent in i?|"^. 

By definition, chirality is reversed by spatial reflection. 
Spatial reflection exchanges a braid in front with its in- 
verse process, a braid behind. It is therefore the same 
as replacing all phases R^^ with (i?"'')*, or equivalently, 
with interchanging Fib^ and Fib*!. Consequently we 
see that any purely diagrammatic model which is based 
on a single Fibonacci UBMTC Fib!|_ and involves non- 
trivial braiding is necessarily chiral (i.e. not reflection- 
invariant). A chirally symmetric model involving non- 
trivial braiding would have to be based on a UBMTC 
involving at least two copies of the fusion algebra Fib. 
An example would be a model based on the UBMTC 
Fib^ X Fib^, which transforms into itself under spatial 
reflection, where the model has a Hamiltonian invariant 
under the exchange of the + and — sectors. 

Note that the chiral nature of the golden chain may be 
concealed by the use of a modified translation operator 
which eliminates all non-trivial braids, for example 



T 



M 





FIG. 25. Structures of operators used in the definition of 
^A,p,D^T ^ j-^j. g^ ^j^g^jj^ q£ p^ qj. ^p^ interacting Fi- 
bonacci anyons. Arrows are not required for Fibonacci anyons 
as all charges are self-dual. 



i ranging from 1 to p inclusive, and identifying i = p + 1 
with i ~ I, but on the disc the constraint applies only 
for 1 < i < p. 

If this restriction is enforced by including terms in 
lT'f_i j j_|_]^ which apply an arbitrarily large energy penalty 

to invalid states, then the behaviour of T^ on the disc is 
such that this restriction is appropriately applied to all 
pairs {i, i-l-l}, 1 < i < p and not to the pair {p, 1}. Thus 
the structure of T° makes it possible to easily enforce the 
restriction on the Hilbert space via Eqs. [12] and [Ml just 
by modifying h^_i j ^,i to impose large energy penalties 
on the unphysical states. 



on the disc. [Applying this operator to the fusion tree 
of Fig. [ijiv), one readily sees that its action may be 
evaluated without requiring knowledge of i?"''.] How- 
ever, without physical motivation for the loop in Eq. IA2I 
this apparent restoration of the chiral symmetry of the 
ground state is artificial, and may be compared to arbi- 
trarily introducing extra factors of —1 into the transla- 
tion operator on a chain of fermions, purely to suppress 
any minus signs (and non-zero momenta) arising from 
fermionic statistics. As a mathematical transformation 
it may sometimes be useful, but the model no longer cor- 
responds directly to the original physical system. 

To understand why the operators T^ and T"^ 
[Fig. [71[iv), Eq. [3 and Fig.[TH] represent the natural def- 
initions of translation on the disc and torus respectively, 
the reader is directed to SecHIl where these operators are 
derived, and also to their application to the Heisenberg 
spin chain in Sec. IIIIFI 



Appendix B: Restricting the Hilbert space 

To restrict the Hilbert space of the spin chain so that 
it yields a valid fusion tree, we must exclude states for- 
bidden by the fusion rules. For example, for a chain 
of Fibonacci anyons we prohibit all states in which Xi = 
Xi+i = 1. On the torus this condition is applied for all x^. 



Appendix C: Construction of the five-body 
operators in Hamiltonian H ' ' ^ 



This appendix presents the construction of the terms 
^n-2,n-i,n.i,2 and /i^_i „.i 2.3 i^ Hamiltonian i/A,p.D->T 



for a chain of AFM or FM interacting Fibonacci 
anyons. These operators will be defined using the fusion 
trees of Fig. (251 1) and (ii) respectively. 

For h^_2 n-i n 1 2 i'^ flis basis of the first diagram of 
Fig- [SHi), we may identify 6i = Xp-i, 62 = Xp, 63 = xi, 
64 = zi , where we have labelled by zi the total charge of 
the five anyons on sites {n — 2, n — 1, n, 1, 2}. The coef- 

/, \ 616263 

ficients I ^„-2.n-i.n,i.2 ) of this operator therefore 

V ■ '/ 6^6^6364 

describe a mapping from states labelled by Xp-i, Xp, xi, 
and zi into states labelled by x'p_i, x'p, x'l, and zi. To 

reproduce the behaviour of the operator hz^_i „ i corre- 
sponding to h^i on the p-site torus, we want this operator 

to be independent of zi, and to correspond to h^_i „ 1 on 

■ ■' ' Its 



the subspace labelled by {xp_i,a;p,a;i,a::' j^, 
coefficients are therefore given by 









V zi 



n — 2,n— l,n.l,2 



(CI) 



21 



Similarly, h^_^ n i 2 3 '^^ ^^^ basis of the third diagram 
of Fig. [SST ii) takes its coefficients from hf, 12 for all 22, 
where Z2 represents the total charge of the five anyons 
on sites {n — 1, n, 1, 2, 3} and we make the identifications 

&1 = Xp, &2 = Xl, 63 = X2, bi = Z2. 

Finally, on the torus the Hilbert space is restricted to 



states which satisfy the constraints ij+i € Xi x t for all 
values of i, and site p + 1 is identified with site 1. On 
the disc this restriction is enforced by the fusion rules for 
1 < I < p — 1 , but for i — p it must be applied manually 
by inserting an arbitrarily large energy penalty into the 5- 
site Hamiltonian terms for states satisfying Xp — xi — 1. 
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